Chaotic Scattering in the Regime of Weakly Overlapping Resonances 
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We measure the transmission and reflection amplitudes of microwaves in a resonator coupled to 
two antennas at room temperature in the regime of weakly overlapping resonances and in a frequency 
range of 3 to 16 GHz. Below 10.1 GHz the resonator simulates a chaotic quantum system. The 
distribution of the elements of the scattering matrix S is not Gaussian. The Fourier coefficients of 
5* are used for a best fit of the autocorrelation function of 5 to a theoretical expression based on 
random-matrix theory. We flnd very good agreement below but not above 10.1 GHz. 
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Chaotic quantum scattering occurs when Schrodinger 
waves are scattered by a system with chaotic classical 
dynamics. For time-reversal invariant chaotic systems, 
the spectral fluctuations of the eigenvalues coincide [l| 
with the predictions of the Gaussian Orthogonal En- 
semble (GOE) of real and symmetric random matrices. 
The eigenvalues manifest themselves as resonances with 
average spacing D and average width F. The theory 
of chaotic scattering has been largely developed in the 
framework of nuclear reaction theory Predictions 
of the theory have been thoroughly tested both in the 
regime of isolated resonances (F <C i?) 0] and in the 
Ericson regime (F ^ D) 0], especially in the context 
of nuclear physics Q but also in several other areas of 
physics 0. In contradistinction, we are not aware of any 
thorough investigation of chaotic scattering in the regime 
of weakly overlapping resonances that would comprise all 
complex reflection and transmission elements of the scat- 
tering matrix. In this Letter, we present data in that 
regime and compare these with theoretical predictions. 

Experiment. We use a microwave cavity made of Cop- 
per coupled to two antennas and measure the response 
to an external field as a function of radiofrequency /. 
The microwave cavity has the shape of a tilted stadium 
billiard 0], see the insert of Fig. [T] The dynamics of 
the classical stadium billiard is chaotic. The tilted shape 
was used in order to avoid bouncing-ball orbits between 
parallel walls. The height of the cavity is 14.6 mm. For 
frequencies / < /max = 10.1 GHz, only a single vertical 
mode in the microwave cavity is excited. In that regime, 
the cavity simulates a two-dimensional chaotic quantum 
system and is a microwave billiard [ij . The experiment is 
performed at room temperature, with Ohmic losses at the 
walls of the cavity. A vector network analyzer couples mi- 
crowave power in and out of the resonator via either one 
or both antennas and yields the complex elements Sabif) 
of the symmetric scattering matrix, where a,b — 1,2. 
The range 3 GHz < / < 16 GHz was covered in steps of 
A — 250 kHz in reflection measurements (yielding 5*11 (/) 
and S'22(/)) and of A = 100 kHz in transmission mea- 
surements (yielding Suif))- Fig. [T] gives examples of the 
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FIG. 1: Absolute squares of the scattering matrix elements 
Sab for signal transmission from antenna 2 to 1 (upper panel) 
and reflection at antenna 1 (lower panel) between 9.0 and 
9.5 GHz. On the logarithmic decibel scale —x dB means an 
attenuation of the microwave power by the factor 10^^^°. The 
resonances overlap and create a pattern of fluctuations. In- 
sert: The shape of the two-dimensional microwave resonator 
used in the experiment. The points 1 and 2 indicate the po- 
sitions of the antennas. 



measured transmission and reflection intensities. 

Fig.[2]shows histograms of the distribution of S'-matrix 
elements in two frequency intervals. The distribution 
of Re{S'ii} is strongly peaked near 1, especially for the 
lower interval, and obviously not Gaussian. The distribu- 
tions of ImlS"!!} and of Re{S'i2} deviate from Gaussians 
(solid lines). The distributions of the phases (rightmost 
panels) are peaked. 

We use the data to construct the S'-matrix autocorre- 
lation functions ^^(e) = Sab{f )S:,{f + e) - \ I' 
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for a, 6 = 1,2. The bar denotes an average over a fre- 
quency window. Three examples for Cab(e) are displayed 
as points in the two upper panels of Fig. [3] (data taken be- 
low /max — 10-1 GHz), and in the insert of Fig. [1] (data 
from above /max where the cavity does not simulate a 
two-dimensional microwave billiard). The values of the 
scattering matrix Sabif) are seen to be correlated, with a 
correlation width F w several MHz. With S'^j = Sab — Sab 
we have also determined the "elastic enhancement factor" 

W = ( l^fiP \s^^\2 \ I i^-fl^p as a function of /, both 

from the autocorrelation functions and from the widths 
of the distributions of the imaginary parts of the scat- 
tering matrix (Fig. [21). Both results agree very well and 
yield a smooth decrease of W with / from W « 3.5 ± 0.7 
for 4 < / < 5 GHz ioW ^ 2.0±0.7 for 9 < / < 10 GHz. 
The computation of the enhancement factors based on 



a theoretical expression for the S-matrix autocorrelation 
function introduced below yields the values W=2.8 and 
W=2.2, respectively. Moreover, we have converted the 
scattering functions Sabif) (measured at M equidistant 
frequencies with step width A) into complex Fourier co- 
efficients Sab{t) with < > 0. Instead of the Fourier index 
k we use the discrete time interval t = k/{MA) elapsed 
after excitation of the resonator. The Fourier coefficient 
Sab{0) is proportional to Sabif)- We find that 5'i2(0) « 0. 
Any two complex Fourier coefficients Sabit) of Sabif) are 
uncorrelated random variables 0]. For t > 0, the co- 
efficients Sabit) have an approximately Gaussian distri- 
bution about their (^-dependent) mean value p^. This 
result is unexpected and was neither predicted theoreti- 
cally nor found experimentally before. 
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FIG. 2: From left to right: Histograms for the scaled distributions of the real and imaginary parts of the reflection amplitude 
Sii and the real part and the phase of the transmission amplitude 5*12, for the two frequency intervals 5-6 GHz (upper panels) 
and 9-10 GHz (lower panels). The scaling factors are given in each panel. The solid lines are best fits to Gaussian distributions. 



The Fourier transform Cabit) of Cafc(e) has Fourier co- 
efficients Xt = \Sabit)\'^- In the lower panels of Fig. [3] 
(in Fig. d]) we show data for log^Q ^ab it) versus t for two 
values of {a, 6} (for {a, 6} = {1,2}, respectively). The 
cutoff at t = 800 US in both figures is due to noise. The 
Sabit) being nearly Gaussian, the distribution Piyt) of 
yt = In Xt is expected to have approximately the form 



Piyt) = exp iyt 







(1) 



where rjt = In xt is given by the expectation value of Xt ■ 
The maximum of Piyt) is at yt = rjt, and Piyt) has a 
strong skewness due to the exponential within the argu- 
ment of the exponential, in agreement with the experi- 
mental data. 



Theory. In the regime of weakly overlapping reso- 
nances, the only theory available is due to Verbaarschot, 
Weidenmtiller, and Zirnbauer [ll[ (in the sequel: VWZ). 
These authors model the scattering matrix 5* of a time- 
reversal invariant system in terms of a GOE Hamilto- 
nian matrix of dimension N . In the absence of "direct 
reactions" (i.e., for Si2iE) = 0), the relevant param- 
eters of the theory are the "transmission coefficients" 
Tc = 1 — I Sccif) P which measure the unitarity deficit 
of the average S'-matrix. Given the Tc, the theory uses 
the limit iV ^ oo to predict for all values of T/D the 



3 



1.0 

0.5 
1.0 





' 1 ■ 




1 ■ 










■ 1 1 1 1 1 I 1 1 





2 

1 — 


4 6 

£(MHz) 

« 1 « 1 


8 10 




800 



FIG. 3: Upper panels: Comparison of the autocorrelation 
function Cabi^-) constructed from the data (points) and the 
fit using Eq. ^ (full line), both normalized by the value 
of C(0) as given by Eq. ([2}. Lower panels: Fourier coef- 
ficients Cab{i) of the autocorrelation functions (points) and 
the Fourier transform of the fit of Cab{^) as given by Eq. Q 
to the data (full line). The elements Si2{f) and Sii(/) were 
taken from the frequency window 9-10 GHz. 
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FIG. 4: Same as Fig. |3]but for the scattering function Si2{f) 
taken in the frequency window 12-13 GHz. 



S'-matrix autocorrelation function 

oc 1 

Cofc(e) = i j d\id\2 y"dA/i(A, Ai, A2)Ja6(A, Ai, A2) 


X exp(-i7re(Ai + A2 + 2X)/D) 
(1 - T,A) 



nTTi 



((l + T,Ai)(l + r,A2))i/2 



(2) 



in terms of the ratio e/D. To simulate Ohmic absorp- 
tion by the walls of the cavity, we introduce additional 
fictitious channels [l^ and associated transmission coef- 
ficients Tc with c = 3,4, . . .. These are defined below. 
The product over channels c extends over both, the an- 
tenna channels and the ficticious channels. The function 
Jab(A, Ai,A2) depends on the As and on the transmis- 
sion coefficients Tq, Tb for the open channels. Both the 
integration measure fi{X, Xi, X2) and Jab are given ex- 
plicitly in Ref. 11]. The correlation width F is actually 
determined by Eq. ([2]) but approximately given by the 
"Weisskopf estimate" T « [D/{2tt)]J2cTc- Equation (ED 
comprises what is known theoretically in the regime of 
weakly overlapping resonances. Higher moments of S 
are not known, not to speak of the complete distribution 
of S'-matrix elements. 

Much more is known both for F ^ D and for F -C -D. 
In the Ericson regime, the distribution of S'-matrix el- 
ements is Gaussian; the correlation function Ca6(e) has 
Lorentzian shape, with F given by the Weisskopf esti- 
mate [l3|; the Fourier transform Cab(t) of Caf,(e) (which 
describes the decay in time of the modes in the cavity) 
is exponential in time; for « 1 (strong absorption) 
the distribution of the phases of the Sat is constant. For 
F ^ Z?, on the other hand, the distribution is far from 
Gaussian. (Consider, f.i., the single-channel case. The 
unitarity condition |S(/)| = 1 confines S{f) to the unit 
circle. The phase of S{f) increases by 2tt over the width 
of every resonance and is nearly stationary in between 
resonances.) The regime F « _D interpolates between 
these two extremes and we expect a non-Gaussian dis- 
tribution of S(/). The results in Fig.[2]give experimental 
information on that distribution and confirm our expec- 
tation. With decreasing /, the distributions deviate ever 
more strongly from Gaussians. As for Cab{t), Eq. ^ pre- 
dicts a power-like decay in time, in striking contrast to 
the exponential decay valid for F ^ £). That prediction 
has been discussed and used in Refs. [l^ and experimen- 



tally tested with microwave resonators in Refs. [13, [H. 
However, these papers did not apply any statistical tests 
based upon a goodness-of-fit (GOF) as done below. 

Analysis. We model Ohmic absorption by a large num- 
ber of absorptive channels with very small transmission 
coefficient each pJi] . The product in Eq. ^ over absorp- 
tive channels is then replaced by an exponential function 
of the sum Tabs of the transmission coefficients of these 
channels, and Ca6(e) depends on Ti,T2,rabs and D. The 
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Fourier transform was fitted to the Xf-data, shown in the 
lower parts of Figs. [3] and IH We used Ta = l — \ Saa P for 
a = 1,2 and calculated the mean level spacing D from 
the Weyl formula This left 

Tabs as the only free pa- 
rameter. In order to allow for secular variations of Tabs, 
the data taken between 3 and 16 GHz were analyzed in 
1 GHz intervals with the help of a maximum likelihood 
fit. Wc find that the sum Ti + T2 + Tabs increases from 
0.11 in the interval 3-4 GHz to 1.15 in the interval 9- 
10 GHz. The resulting increase of Tabs is consistent with 
conductance properties of Copper. Using the Weisskopf 
estimate we find that T/D increases from 0.02 to 0.2 over 
the same range. This shows that we deal with weakly 
overlapping resonances. The results of the fits are shown 
as solid lines in the lower two panels of Fig. [3] and in the 
lower panel of Fig. |3) For an exponential decay in time, 
the curves in these panels should be straight lines. This 
is clearly not the case. The solid lines in the upper two 
panels of Fig. [3] and in the upper panel of Fig. |4] are the 
Fourier transforms of the VWZ fits. In Fig. [3] they agree 
well with the data points, up to small discrepancies which 
are attributed to finite-range-of-data errors. In the up- 
per panel of Fig. |4]the discrepancy between fits and data 
points is displayed more clearly than in the lower panel. 

The quality of the agreement between data and fits 
in Figs. [3] and |4] is assessed in terms of a highly sensi- 
tive goodness-of-fit (GOF) test (the Fourier coefficients 
scatter over more than five orders of magnitude!). The 
fit of Eq. ^ determines the expectation value of Xt 
and, thus, r]t = Inxt in Eq. ((!]). If the distribution of 
the ut were Gaussian, the GOF test would be defined in 
terms of J^tiVt — VtY ■ The appropriate generalization 
for the distribution P{yt) in Eq. ([T]) is the expression 
/ oc [exp (yt - ?7t) - (yt - ?7t) - 1], see Chaps. 14, 16 
of Ref. [17| . This quantity is non-negative and vanishes 
exactly if the data coincide with the model for all t. For 
large M, / is approximately x^~distributed with M de- 
grees of freedom. For each frequency interval of length 
1 GHz we have M — 2400, since each of the three ex- 
citation functions 5*11 (/), 5'i2(/) and 5*22 (/) contributes 
800 Fourier coefficients. We admit a 10 % probability 
for an erroneous decision. The fit using Eq. ^ is ac- 
cepted in all intervals below / = 10 GHz and is rejected 
in all intervals but one above 10 GHz. A similarly thor- 
ough and mathematically reliable test of the theory of 
chaotic scattering has not been performed before, see 
Refs. [13, [3]. This fact motivated our work. We con- 
clude that Eq. ([2]) is compatible with our data as long as 
the resonator supports only two-dimensional modes and 
simulates a chaotic billiard. We have numerically simu- 
lated the fluctuations above 10.1 GHz under the assump- 
tion that the two vertical modes do not interact and the 
Hamiltonian matrix is block-diagonal, each block taken 
from the GOE. In this way we reproduced qualitatively 
the results of Fig. 01 In the sense that the GOE describes 
full chaos, a block-diagonal random matrix represents 



additional symmetries. The disagreement between the- 
ory and experiment above 10.1 GHz shows that our test 
is sensitive to the existence of such symmetries. We con- 
clude that first, in the regime of overlapping resonances 
our test is sensitive to symmetries in a Hamiltonian sys- 
tem and second, that Eq. ^ is compatible with the data 
as long as the scattering system is fully chaotic. 

Summary. We have investigated a chaotic microwave 
resonator in the regime of weakly overlapping resonances 
r « Z?. The distributions of S'-matrix elements are not 
Gaussian. In each of 13 frequency intervals we deter- 
mined 2400 uncorrelated Fourier coefhcients of the ele- 
ments of the scattering matrix. Surprisingly, these have 
nearly Gaussian distributions. The data were used to test 
the VWZ theory of chaotic scattering. The predicted 
non-exponential decay in time of resonator modes and 
the frequency dependence of the elastic enhancement fac- 
tor are confirmed. Our goodness-of-fit test is based on 
a large number of data points and constitutes the most 
sensitive test of the theory of quantum chaotic scattering 
for weakly overlapping resonances performed so far. We 
show that VWZ is compatible with the data as long as the 
resonator simulates a fully chaotic quantum system. The 
theory can, thus, be used with confidence to predict av- 
erage cross sections and S'-matrix correlation functions. 
The agreement fails when a second vertical mode ap- 
pears. This suggests that our analysis may serve as a 
tool to detect symmetries and/or regular motion within 
a chaotic system in the regime of overlapping resonances. 
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